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Abstract 



We show that S.Vavasis' sufficient condition for global invertibility of 
a polynomial mapping can be easily generalized to the case of a general 
Lipschitz mapping. 

Keywords: Invertibility conditions, generalized Jacobian, nonsmooth 
^ ' analysis. 

In applications, in particular in finite-element analysis it is useful to have 
some sufficient conditions for testing invertibility or injcctivity of a mapping 
/ : R™ — » R". Usually /(•) is defined on some simple region /, say cube or 
simplex, in the reference domain. The image /(J) is some region with more 
\ complicated geometry (a grid cell) in the physical domain. 

■ In a paper |jj S.Vavasis proposed a sufficient condition for global invertibility 
of a polynomial mapping in R™. Here we show that this condition can be more 
naturally formulated in the framework of nonsmooth analysis (see, e.g. 0), 
and this enables to spread the results of G[ and translate them almost word- 
for-word into a seemingly much more general setting (at least in the linguistic 
sense). Thus this note may be regarded as a feedback to Q]. 

We need several definitions. The standard reference on nonsmooth analysis 
is M . Our discussion is restricted to a finite-dimensional case though an essential 

■ part of H is devoted to the infinite-dimensional generalizations. 

Recall that nonsmooth analysis works with Lipschitz functions that are al- 
most everywhere differentiable. Roughly speaking, the generalized (or Clarke's) 
gradient dF of a Lipschitz function /(•) : R n — > R at the point x £ R n is 
defined as a convex hull of (almost) all converging sequences of the gradients: 

df(x) d = conv(\imV(f(xi)), for Xi — > x and /(•) is differentiable at points Xi 
and the sequence Vf(xi) converges. It is essential that at the points of smooth- 
ness of /(•) the generalized gradient coincides with gradient, and for a convex 
function — with its subgradient. 

Similarly the generalized Jacobian dF of a Lipschitz mapping F(-) : R n — > 
R m at the point x is a convex hull of all m x n matrices obtained as limits of 
sequences JF(xi) (its Jacobian matrix at the point of smoothness Xi), where 
Xi — > x and F(-) is differentiable at a;,-. 

'Supported in part by RFBR grant 02-01-00716 

1 Actually, Clarke's definitions of the generalized gradient and the generalized Jacobian are 
somewhat less restrictive. We may assume that the points of the sequences Xj or JF(xi) 
involved do not belong not only to nonsmooth point of the map but additionally they do not 
belong to an arbitrary set of Lebesque measure zero. Such definitions are technically more 
convenient. 
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The generalized Jacobian has maximum rank at point xq if each matrix from 
dF(xo) has maximum rank. 

Nonsmooth inverse function theorem |^|, Th.7.1.1] states that Lipschitz map- 
ping F : R n — > R n , whose generalized Jacobian has maximum rank at point xo 
is locally Lipschitz invertible. 

In [Q S.Vavasis proposed a sufficient condition for the nondegeneracy of a 
matrix family A4 consisting of square n x n matrices. To present this result in 
a slightly more general form we need some (more or less standard) definitions. 

A cone K, C R n is a set with a property that for all x £ K, and for any 
A > 0, Aa; G /C, i.e. K, contains all intersecting rays through the origin. 

A cone K. is convex if sum of any points of K, belongs to KL. 

A convex cone is polyhedral it can be represented in the form: 1C — {x\x — 
E7=i A ^ > 0, a 4 G R™}. 

A convex cone K, is acute if it contains no nontrivial subspaces or equivalently, 
if no finite set of elements of K. sums to zero. 

If A C R N then cone hull (A) is a union of all rays through the origin 
intersecting A. 

Let formulate Vavasis' condition. Denote by Ri C R™ the set of all i-th 
columns of the matrices belonging to M. and assume that all Ri , i = 1, . . . , n 
are separated from the origin, i.e. for some 8 > Ri fl {x | ||x|| < 8} = (here 
|| ■ || denotes the Euclidean norm). Set Ki = cone hull(Ri), i = 1, . . . , n and let 
any cone out of the 2 n_1 cones K\ ± K% ± . . . ± K n be acute. 

By definition, V -family is any matrix family satisfying these conditions. It 
follows from the above that all matrices in ^-family are nondegenerate. 

If additionally, any cone involved K\ ± Ki ± . . . ± K n is contained in a cone 
K a = {x G R™ | ax > e\\x\\} for some certificate vector a G R™ and e > 
then such family is defined as strict V -family. In particular, if all Ki (or Ri), 
i = l,...,n are polyhedral then by Farkas lemma all V- families are strict. 
Moreover, checking that some matrix family is V^-family is reduced to solving 
2™ _1 linear programs, and thus the overall test could be performed in linear 
time with respect to the input in any fixed dimension. 

Theorem 1 Let F(-) = (/ J (-), /"(•)) : U G R™ — > R" be any Lipschitz 
mapping defined on the convex reference domain U . If the set of the generalized 
Jacobians J = { J G dF(u), u G U} forms a strict V -family then F(-) is globally 
invertible on U . 

Proof. Actually the demonstration is a direct translation into the new 
setting of the original proof from Q| . 

At first, local invertibility of F(-) at any point u £ U of the reference domain 
follows from the nonsmooth inverse function theorem as the generalized Jacobian 
dF(u) has maximum rank at any point by construction. 

Secondly, to show global invertibility it is enough to check injectivity of 
the mapping. Take any different points u, v G U from the reference domain. 
By convexity of U and by the nonsmooth analog of the Lagrange formula 
Th.2.6.5] for almost all pairs u,v G U the following equality holds: F(v) — 
F(u) — J Q JF(u + t(v — u))(v — u)dt. Assume w.l.o.g. that the first coordinate 
of the vector v — u is nonnegative (otherwise, exchange v and u). Set Ki — 
cone hull(df l (x), x G U). Now assume for a moment that v — u > 0. Take the 
corresponding certificate vector a for the cone K\ + . . .+K n . By construction, for 
all s eU, (a, dfi(s)) > eS, i = l,...,n. Hence, (a, F(v)-F(u)) > eS\\u—v\\ > 
and the injectivity follows. In general case, we take any certificate vector for 
the cone K\ ± K2 ± . . . ± K n with the same sign pattern as the sign pattern of 
the coordinates of the vector v — u. 
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Informally, Vavasis sufficiency condition assumes that the columns of the 
corresponding Jacobian matrices are independent and thus may be extremely re- 
strictive but, on the other hand, it has attractive enough decomposition feature. 
Namely, let the reference domain I be subdivided into several simple regions, 
say, I is a cube that is partitioned into parallelepipedal patches / = Uf/^P by 
several axis- parallel hyperplancs. Now let P(-) = (f 1 ^), ■ ■ . ,/"(•)) : I ~ * R n 
be any Lipschitz mapping. Let Aj = \J xe i i JF(x) (the union of all generalized 
Jacobians in the closure of the patch ij). Obviously these sets can be computed 
separately for any closed patch /j i = 1, . . . , M. As above we obtain the following 
sufficient condition: if each of the sets A,-, i = 1, . . . , M forms a strict V- family 
and there exists a certificate vector common to all Aj then the mapping P(-) 
is globally invertible. Equivalently, as Vavasis condition is insensible to taking 
convex hulls, F(-) is globally invertible on / if the set A = conw(U*£ 1 Aj) forms 
a strict V-family. 

For a simple example, assume that P(-) is a continuous piecewise polynomial 
(product) Bernstein-Bezier (BB) mapping as proposed in , i.e. in any patch 
Iq, q = 1, . . . , M the mapping is given by 

f&, ...,£«) = £•••£ n u ... in ( p )^ (i - . . . - e«rs 

j 1= i n =0 ^ X ' 

where (£i, . . . , £„) are local coordinates in I q = {0 < £i < 1, . . . , < £„ < 1} 
and the vectors /,? i £ R", ii, . . . , i n — 1, . . . , n are called the control points. 
It is well known that the image F(I q ) of BB mapping is contained in the convex 
hull of the control points. It is also well known that the partial derivatives of 
F(-) can be put into BB form (and the resulting control points for the derivatives 
are effectively computable linear combinations of the control points for F(-)). 
Thus for any i = 1, . . . ,n the set of all i-th columns of the Jacobian matrices 
JF(I q ) — {JF(x), x S I q } is contained in some explicitly computable polytope 
Pf . Hence, all possible i-th columns of the Jacobian matrices JF(x), x E I 
are contained in Pj = conv{P} , . . . , Pj M }, i — 1, . . . , n. (Here we use the 
continuity of F(-) as at any point belonging to the intersection of some patches 
x € lj! H . . . n Ii k the generalized Jacobian JF{x) is by construction contained 
in conv{JF(Ii 1 , . . . , JF(Ii k )}.) Thus if the matrix family V, whose z-th column 
belongs to Pj, i = 1, . . . , n, form a strict V^-family then P(-) is globally invertible 
on the whole domain I. 

Here we show that invertibility test in |jj is valid for a much more ample class 
of mappings and thus should be simultaneously robust and restrictive enough. 
Therefore it would be nice if some new arguments would be applied to the 
invertibility problem even in the simplest case of bilinear polynomials. 
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2 The reader should take into account that the example below is almost explicit in |lj 
and could be easily recovered from the arguments therein but with some additional formal 
arguments. 
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